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Abstract. Wc address the question of whether numerical particle-in-cell (PIC) 
simulations and laboratory laser-plasma experiments can (or will be able to, in the 
near future) model realistic gamma-ray burst (GRB) shocks. For this, we compare 
the radiative cooling time, t coo i, of relativistic electrons in the shock magnetic 
fields to the microscopic dynamical time of collisionless relativistic shocks — the 
inverse plasma frequency of protons, uipp . We obtain that for t coo \u)pp < few 
hundred, the electrons cool efficiently at or near the shock jump and are capable 
of emitiing away a large fraction of the shock energy. Such shocks are well-resolved 
in existing PIC simulations; therefore, the microscopic structure can be studied 
in detail. Since most of the emission in such shocks would be coming from the 
vicinity of the shock, the spectral power of the emitted radiation can be directly 
obtained from finite-length simulations and compared with observational data. 
Such radiative shocks correspond to the internal baryon-dominated GRB shocks 
for the conventional range of ejecta parameters. Fermi acceleration of electrons in 
such shocks is limited by electron cooling, hence the emitted spectrum should be 
lacking a non-thermal tail, whereas its peak likely falls in the multi-MeV range. 
Incidentally, the conditions in internal shocks are almost identical to those in laser- 
produced plasmas; thus, such GRB-like plasmas can be created and studied in 
laboratory experiments using the presently available Petawatt-scale laser facilities. 
An analysis of the external shocks shows that only the highly relativistic shocks, 
corresponding to the extremely early afterglow phase, can have efficient electron 
cooling in the shock transition. We emphasize the importance of radiative PIC 
simulations for further studies. 



PACS numbers: 98.70.Rz, 52.72,+v, 52.35.-g, 95.30.Qd, 52.65.Rr 



Submitted to: Astrophys. J. 



1. Introduction 



It has been shown in recent years that collisionless relativistic shocks are mediated by 
the Weibel instability — a current filamentation instability that produces strong, sub- 
equipartition magnetic fields at the shock front (Mcdvcdcv & Locb 19991. This is an 
attractive model for gamma-ray bursts (GRBs), because it puts a synchrotron shock 
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model on a firm physical ground. However, it has been thought that the generated 
magnetic fields should occupy an extremely small volume of the shocked region, hence 
the radiation emitted by the relativistic electrons in these fields should be extremely 
weak and astrophysically unimportant. Conventionally, one assumes that the strong 
Weibel fields live on scales of tens of plasma skin depths (tens of c/lo pp ) behind the 
shock, which ranges from only few centimeters in internal shocks to about 10 cm 
or so in external shocks; both are some ten orders of magnitude smaller than the 
typical size of the shocked region. Therefore, unless the synchrotron cooling length, 
^cooic, is comparably small, the radiative efficiency of such shocks should indeed be 
very small. This further raises a concern that a number of numerical particlc-in- 
cell (PIC) simulations attempting to model GRB shock^j] (Spitkovsky 2008 Chang 
Keshet et al 2008) probe spatial scales that are too small, and, hence, 



et al 2008 



cannot be used to deduce the radiation spectrum of GRB emission. This also can 
cast some doubt on jitter radiation as a viable explanation of spectral properties of 
prompt GRBs (|Medvedev||2000[ |2006| [Hededal||2005| ). R ecently, long dura t ion shock 



al 



2008 



simulations (Keshet et al 20081 and analytical works (Katz et al 



ISironi fc Goodman||20071 |Goodman fc MacFadyen[ [2007 



2007 Couch et 



Milosavljevic et 



al 20071 have suggested that appreciable magnetic fields could still survive in the 
downstream region due to the effects of accelerated particles or the field amplification 
by upstream turbulence. Regardless of the ultimate resolution of the magnetic decay 
question, there exists another physical regime where most of the radiation is produced 
near the shock, rather than in its downstream. Simulations of such shocks then capture 
most of the relevant physics without the need to simulate the far downstream region. 
In this paper, we explore the relevance of such shocks to GRB scenarios. 

We use the results of recent numerical simulations of relativistic collisionless 



shocks in electron-ion (Spitkovsky 20081 and electron-positron (Chang et al 2008 



Keshet et al |2 008 ) plasmas to estimate and compare the electron cooling and plasma 



time-scales in the shock transition region. We then deduce the parameter range of 
internal (both baryon-dominated and electron-positron pair-dominated) and external 
shocks, for which the shocks can be considered as radiatively efficient, i.e., the radiated 
energy in the shock transition is of order the electron energy j e m e c 2 , where 7 e is an 
average Lorentz factor of the bulk electrons. We find that internal shocks are in this 
radiative regime for a reasonable range of model parameters, V ~ 200 — 300 and 
Tj > 2.5, where V and I\ are, respectively, the bulk Lorentz factor of the ejecta and 
the Lorentz factor of an internal shock inside the ejecta, measured in the center of 
mass frame of the colliding shells. In contrast, external shocks may be in the radiative 
regime only at very early times, less than 100 seconds after the explosion. We also 
argue that the physical conditions in internal shocks are very close to those that can be 
achieved in laser-plasma interaction experiments. Therefore, the physics of the Weibel 
instability and jitter radiation can be studied experimentally at laser-plasma facilities, 
such as Omega EP, NIF and others. In particular, the Hercules experiment has been 
proposed and is now under development at the University of Michigan ( Reynolds et al 
|2006||2007| |Huntington et al |2008"| for technical details, see |Maksimchuk et al |2008[ ) 



to create and diagnose the Weibel instability and turbulence in the laboratory high- 
energy density plasmas, as a part of the Laboratory Astrophysics and High-Energy 
Density Physics programs. 



§ These 2D simulations studying shocks in electron-position plasma are extending to at most ten 
thousand skin depths, i.e., of order tens of meters (in the shock co-moving frame) for internal shocks. 
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2. Dimensionless electron cooling time 

The electron synchrotron cooling time is obtained from the relation that the emitted 
energy is of order the initial electron's energy: 

-Psyntcooi = (4/3)o- T c7e (B 2 /8ir)t cool ~ 7 e m e c 2 . (1) 

Here, <tt is the Thompson cross-section, 7 e is the Lorentz factor of the emitting 
electron, and P syn is the synchrotron emission power per electron. This gives 



Note that this expression holds for jitter radiation as well, because the total emitted 



power by an electron in the jitter regime is identical to that of synchrotron ( Medvedev 



2000) 



We will compare t coo \ to the characteristic time of the plasma processes at a shock, 
the inverse relativistic plasma frequency of the protons, u>~ p l , which is 

^ , (3) 

where n' is the particle density behind the internal shock measured in the downstream 
frame. To within a factor of two, the so-defined oj pp corresponds to the nonrelativistic 
upstream plasma frequency, because n' = 4nPi, where n is the particle density of 
the unshocked ejecta in its own co-moving frame. It is natural to introduce the 
dimensionless cooling time 

-Fcool — ^cool^pp- (4) 

In our analysis we will use the results of electron- ion (Spitkovsky 2008) and 



electron-positron (Chang et al 2008) shock PIC simulations. Apparently, the 
evolution of currents and magnetic fields in the vicinity of a shock in electron-proton 
and electron-positron plasmas does not differ substantially, because in the electron- 
proton case the electrons in the downstream carry about 30% — 50% of the proton 
energy and, hence, their effective mass is comparable to the proton massf]]] 

Figure [T] represents a snapshot of a steady state shock in the electron-ion 
plasma obtained from a PIC simulation with Lorentz factor T = 15 and mass ratio 



mi/m e = 100 (|Spitkovsky 2008[ ). This and other simulations (Chang et al 2008 



Keshet et al 2008) indicate that the strongest magnetic fields with energy density 
> 10% of the kinetic energy density occupy the region of few tens of ion plasma 
skin-depths around the main shock compression, as measured in the frame of the 
downstream fluid. Hereafter, we will use the value of I ~ 50c/u> pp for the transition 
length (for electron-positron shocks, the ion plasma frequency is replaced with the 
plasma frequency of pair plasma). Since the fields are highly inhomogeneous, we 
assume a factor of two uncertainly in this number. The shock moves at v = c/3 (or 
c/2 in 2D) in the downstream frame, hence the electron residence time in the region 
of high field is t rcs ~ l/v ~ lhQuj~ p l . Taking the uncertainty in / and other parameters 
into account, we estimate that t ICS ~ (100 — 300)0^. We did not account here for 

| One has to be careful here because the parallel and perpendicular relativistic masses are 
not identical, m et ± ~ ■Je'me ~ I^mp and m e it ~ ~ Tim p (Tim p /m e ) 2 2> Tim p . 

Hence, the electrons will be effectively much heavier than the protons with respect to their 
acceleration/deceleration, but should behave similar to the protons with respect to deflection and 
pitch-angle scattering. 



Radiative cooling in relativistic shocks 



4 



the dumpiness of magnetic inhomogencities, which shortens the effective i res , and the 
electron trapping in high-field clumps, which is increasing the effective residence time. 
Interestingly, in the far downstream region, the magnetic energy decreases with the 



distance d from the shock as oc d (Chang et al 2008). An electron in this region 



loses energy logarithmically slowly, E oc ln(d). If this scaling holds through a large 
distance downstream, this implies two possibilities. First, if i coo i < t ros , then the 
electrons lose their energy quickly near the shock jump; hence the radiative efficiency 
of such a shock is high. The electron cooling will substantially affect the structure of 
the downstream region, in this regime. Second, if t coo \ 3> t TCS , then the electrons lose 
only a small fraction of their energy, which makes this shock radiativcly inefficient. 
Therefore, we refer to a shock as a "radiative shock" if T coo i < 300 and as a "weakly 
radiative shock" otherwise. Interestingly, if T coo i < f — 50, radiative cooling might 
be substantial already in the upstream region (where the Weibel instability creates 
current filaments and magnetic fields), before the main shock compression. We expect 
that the very formation of such a shock may be strongly affected by radiative cooling 
of electrons. We refer to the regime of T coo i < f — 50 as the "strong cooling regime." 

3. Shock models 

3.1. Internal shocks in bary on- dominated ejecta 

Internal shocks occur inside the ejected material when inhomogeneities, often referred 
to as "shells' that have different masses and velocities, collide with each other 
relativistically. Radiation emitted at the internal shocks is thought to be the prompt 
gamma-ray radiation of a GRB (e.g., Meszaros||2006[ ). 



Here we conventionally assume that a central engine produces an ultra-relativistic 
wind with the kinetic luminosity, L ~ 10 52 erg/s and the Lorentz factor T > 100 
flMeszaros fc Rees||2000| |Rees fc Meszaros||2005| |Pe'er fc Waxman||2004[ ) . At a radial 



distance, R, from the central engine, the co-moving particle density of the ejecta is 

n = 47rR2 ^ mpc3 - (1.8 x 10- cm-3) i52 ^r-, (5) 

where we use the convention that L52 = L/(10 52 erg/s) and similarly for other 
quantities, whereas if no numerical subscript is present, then the quantity is in CGS 
units. The density in the downstream of an internal shock is 

ri = 4I>. (6) 

The magnetic field in the downstream and the electron bulk Lorentz factor are 
calculated as fractions cb ami e e of the post-shock thermal energy density, U = 
n' '('ypmpC 2 + 7 e m e c 2 ) ~ n'TiiripC 2 (in the last equality, we neglected the electron 
contribution, for simplicity, though it may change U by a factor e e /[l + e e ]), 

B' = (87rr 4 m p c 2 n'e B ) 1/2 ^ (1.6 x 10 7 G) L^T^ 1 e X J 2 . (7) 

7e = (m p /m e )T t e e ~ 1.8 x 10 3 rV e . (8) 

For the parameters e e and cb one has to use the emission- weighed quantities, 
rather than the conventionally used volume-averaged ones, because fields are highly 
inhomogeneous in the shock region. This effect is particularly important for es and 
less so for e e because high-energy electrons are distributed more or less uniformly. PIC 
simulations indicate that e e ~ 0.5 (electrons carry 50% of the proton energy) and that 
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e_e reaches unity locally ( Spitkovsky||2008 Chang et al [2008 1. High energy radiation 



from electrons is not included in the present code. In order to calculate the emission- 
weighted e e and £b (black curves in panels (d) and (e) in figure [TJ hereafter denoted 
with a "hat" symbol) we took into account that radiative losses are proportional to 
B 2 j 2 . Since the electron cooling occurs predominantly in the high-field regions, a 
larger effective £b and e e would be deduced from observations. Taking into account 
the factor of 8 difference in the definition of es here and in Spitkovsky 2008; Chang 
et al |[2008 we estimate the emission-weighted value as cb ~ 1/8 ~ 0.1. 



Finally, the dimcnsionlcss cooling time in baryon-dominated internal shocks 
becomes 

TiV" ~ mL-^T 2 Tf Rl2 e- B H-\ (9) 

We assume that the shells have different masses but carry similar linear momenta, 
p ~ T s m s c ~ T r m r c. That is, in this scenario, a constant driving force is assumed 
to continuously eject material while the central engine is active. Thus, the variation 
in the Lorentz factor of the shells is simply reflecting the variation of the density of 
the ejecta (i.e., that of the shell masses). An attractive feature of this model is that 
the total energy of the ejecta is uniformly distributed among individual shells, which 
nicely fits into the relativistic wind picture, unlike the model with shells of equal mass, 
in which the fastest shell essentially dominates the energetics of the outflow. In the 
collision of two shells, two shock waves are formed, which propagate through each 
shell ( Kobayashi et q?"||1997 ). The center of mass Lorentz factor of the two shells is 



„ m,r, + m r T r 

r c = — — (io) 

[m 2 + m 2 + m s m r (T s /T r + r r /r s )] 

so r c ~ y2T s for T s m s ~ T r m r and T r 3> T s . The thermal energy density in the 
downstream of each shock is related to the Lorentz factor of the shells in their center 
of mass frame: 

(Ti) a , P ~ (i/2)(r s , r /r c + r c /r s , r ), (11) 

so that the Lorentz factor of an internal shock can be as high as T { ~ 2" 3 / 2 r r /r s , 
which is at least a few if the Lorentz factors of the "slow" and "rapid" shells are, say, 
T s ~ 100 and T r ~ 1000, respectively, that is m s /m r ~ 10. 

We assume that the central engine is an accretion disk around a solar mass black 
hole. Thus, the flow forms at a few to ten Schwarzschild radii. This sets the variability 
time scale, t v ~ 10~ 4 s. Present observations confirm the variability on a millisecond 
scale, but arc not capable of resolving sub-millisecond time scales. In this model, the 
relativistic shells are ejected from the radius Ro ~ ct v and their initial separation is 
of order i? as well. A collision of the two shells will occur at the radius 

Ri ~ 2R (T- 2 - T" 2 )" 1 - 2T 2 s ct v - (6 x 10 10 cm)L^,,_ 4 - (12) 

One should keep in mind that collisions of the shells can occur at substantially smaller 
radii, if the wind modulation occurs due to instabilities in the outflow (e.g., a jet), at 
radii R^> Rq where the flow is already ultra-relativistic ( Rees & Mcszaros 2005 ) . 



Due to the presence of electrons in the baryonic outflow, the co-moving optical 
depth due to Thompson scattering, Tj, ~ hctt^R/V), approaches unity at the radius 

R ph = T/(na T ) ^ (1.2 x 10 13 cm)i 52 L 2 : 3 , (13) 

where ax is the Thompson cross-section and R/T is the co- moving length. For a 
radiative luminosity of a GRB of L 7 ~ 0.1L (we assumed a 10% radiation efficiency), 
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the electron-positron pair opacity is obtained from the balance between the rates 
of annihilation and of pair production. The latter is the photon density above the 



electron rest mass over the co-moving dynamical time (Pe'er & Waxman 2004). 
Taking into account that only the photons with energies above m e c 2 can produce 
pairs, the pair-producing luminosity is L> ~ (fi,i/ syn /m e c 2 ) /3 L 7 if hv ayn < m e c 2 
and L> = L 7 otherwise, where v syn = j 2 (eB' /m e c) is the co-moving synchrotron 
frequency, (3 ~ 0.5 is the high-energy spectral slope vF v oc . Finally, we obtain 
the pair optical depth to be equal to the square root of the co- moving compactness, 
t± ~ Z' 1 / 2 ~ (L 7 "<TT/47ri?r 3 m p c 3 ) 1 ^ 2 . The radius of the pair photosphere, t± ~ 1, is 

i?pair = i>CT T /47rr 3 m p c 3 ~ (2.2 x 10 3 cm)L 7i51 r^ 3 . (14) 
In the last expression we omitted the spectral correction, for simplicity. Note that non- 



thermal spectra can be seen even for relatively large t± of few tens (Pe'er & Waxman 



2004 1 . Since the optical depth increases with decreasing radius as t± oc l njl oc R 1/2 



one can observe non-thermal spectra produced at internal shocks at radii as small as 
~ 10~ 3 i?p a i r or so. Note also that the position of the pair photosphere depends on 
(unknown) radiative efficiency. 

Figure [2] shows the regions of strong and weak radiative cooling versus R and 
r for a reasonable and a rather extreme cases of Ti = 4 and Ti — 10. One can 
see that internal shocks in the strong cooling regime , T coo i < 10, can occur for 
T < 10 2 1 - 125 and at R < few x (10 10 - 10 11 ), well below the photosphere at such 
low T's; hence, radiation from such shocks is hardly observable. Strongly radiative 
shocks T coo i ~ 100 — 300 can occur above the baryonic photosphere for reasonable 
values of r < 200 — 300 and Ti > 2.5; hence they are likely observable. At the 
photosphere, i? p h = Ri, we have 

R hph ~ (4.9 x 10 11 cm)L%hl% (15) 

T(R hph ) ~ 290i^ 5 t" 1 _ / 4 5 ) (16) 

which depend weakly on the central engine parameters. Excluding R and T in Eq. 
Q9J) , we obtain 

TL P \R^) - 250Ll{ w tl%T 7 %^-\ (17) 
We plot the the regions of strong and weak cooling versus R and in the left panel 
of figure [3] and r coo i vs. Ti in right panel. We conclude that internal shocks with Ti 
as low as 2.5 can be radiatively efficient. 

The peak of the synchrotron radiation in the observers frame is 

u syn = (eB/m e c)j 2 e T ~ 400 MeV L^R^T? e^ /2 e 2 . (18) 
If radiation is emitted in the jitter regime, the spectral peak is expected to be a bi t 



higher: Vj ~ v syn \j '{mi/m e )eB ~ Wv syn (Medvedev 2000 Medvedev et al 2007) 



Thus, for the assumed parameters the emission is expected to be peaked at few tens 
to few hundreds MeV, if the shock is in the radiatively efficient regime with T ~ 300 
and R - 10 12 cm. 

3.2. Internal shocks in pair- dominated outflows 

This case can be readily obtained from the baryon-dominated case by replacing m p 
with m e . We obtain 

T^oti - 5-9 x 10 & L- 2 1/2 T 2 T- 3 R 12 e B 1 i- 1 . (19) 
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Thus, internal shocks in the e -dominated ejecta have a very long radiative cooling 
time and, hence, are radiatively inefficient in our language. The peak of the 
synchrotron radiation in the observers frame is 

v syn ~ 0.12 keV Ll^R^Tp^el (20) 
3. 3. External shocks of afterglows 

In the post-shock region of an external shock propagating into an external medium of 
density n cxt , we have 

n' = 4Tn cxt = (400 cm- 3 )r 2 n cxt , (21) 

7e = (m p /m e )Te e = 1.8 x 10 5 r 2 e e , (22) 

B' = {^Tm p c 2 n'e B ) 1/2 ~ (39 G)T 2 n 1 c i 2 e 1 J 2 , (23) 
Using these quantities, we estimate the cooling time as 

T^S - 7-3 x 10 3 r 2 - 3 rCxt /2 ^e (24) 
This relation is plotted in figure [|for a very dense external medium, n cxt ~ 100 cm 3 . 
Even for such an extreme case, the external shock is radiative only for very large 
Lorentz factors T > 500. It is an order of magnitude larger for a more conventional 
value of rt cxt ~ 1 cm -3 . 

We now consider two models of the external medium density profile, namely the 
constant density interstellar medium (ISM) and the wind outflow models. 

In the ISM model, the Lorentz factor of a blast wave depends on the observed 
time as ( Granot et al ||1999 l 



r^350E 5 1 /Vs^(l + ^) 3/8 *- 3/8 , (25) 
where E ~ ££qrb is the isotropic energy equivalent of the blast wave, <grb is the 
duration of a GRB and n cxt = hism = const is the external medium density. In this 
model, the cooling time is 

T AG,ISM ^ 180 £- 3 3/8 n -V8 (1 + ^-9/8^1^9/8. (2g) 

In the Wind model, the blast wave is propagating in the wind environment with 
the density decreasing with distance as n cx R~ 2 \ the Lorentz factor of the blast wave 
and wind density are ( Chevalier fc Li|[2000 1 



r ~27QEl /2 A- 1/A {\ + zf/H- 1 '*, (27) 

n cxt ~ (1.1 x 10 5 cm- 3 )^ 1 ^! + z)*" 1 , (28) 

where A* = [M W /(1Q- 5 M yr~ 1 )]/[VV/(10 3 km s -1 )] is the wind parameter, M w is 
the mass loss rate, Vw is the wind velocity. The cooling time is 

T C A «< W c 1.2E^A-^(1 + z)S'H- B H-W\ (29) 

In figure [5] we plot T coo i in the external shock versus time after the burst for 
both models and for two sets of afterglow parameters. We see that, except for the 
very early times, in both models, the emission from external shocks should be coming 
from far downstream, not from the main shock compression region. However, in the 
Wind model, the external shock can be radiative up to ~ 100 s after the burst, whereas 
for the ISM model, the radiative shock regime can occur only at times earlier than a 
second after the explosion. Since the afterglow usually sets in at least several tens of 
seconds after the explosion, we conclude that only very early afterglow emission can, 
in principle, come from radiative shocks and only in the Wind model. 
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4. Discussion 

We compared the electron cooling time to the dynamical microscopic time (the plasma 
time, o-ip 1 ) for the internal shocks in both electron-positron and baryon-dominated 
relativistic GRB outflows and for the external afterglow shock. We used the most 
recent PIC simulations to obtain e e , cb and the size of the region with strong magnetic 
field. We evaluated and used in our analysis the emission-weighted e_g, instead of the 
conventional volume-averaged quantity (for e e both emission- and volume- averaged 
quantities are expected to be similar). We then evaluated the residence time of 
an electron in the high-field region as the size of this region over the flow speed. 
Thus, we did not account for dumpiness and less-then-unity filling factor of magnetic 
inhomogeneities (which is decreasing the effective residence time) and the electron 
trapping in high- field clumps (which is increasing the effective residence time). Using 
PIC simulations, one can study the radiative cooling effects much more accurately. 

We obtained that if T coo i = t coo \Lo p < (100 — 300), an electron has enough time to 
radiate away energy comparable to its initial energy, ~ 7 e m e c 2 . If (10 — 50) < T coo i < 
(100 — 300), much radiation is emitted from the strong field region at the shock jump. 
Therefore, we call this regime the "radiative shock" regime. If T coo i ;$ (10 — 50), 
the electron cooling is extremely fast and strong radiative losses are expected in the 
upstream region as well; hence we refer to it as the "strong cooling regime" . Note 
that in both these cases, Fermi acceleration of electrons is hardly possible because the 
electrons lose their energy every time they cross the shock. Therefore, we expect that 
the electron distribution will not be a power-law and, hence, a hard non-thermal tail in 
the radiation spectrum is hardly produced. If T coo i 3> 300, radiation from the shock 
(if any) will be coming from the far downstream region. This region has not been 
fully analyzed in simulations, so it is too early to draw any firm conclusions about this 
regime. However, if we extrapolate the decay rate es oc l/d {d is the distance from 
the shock front) as seen in PIC e ± simulations, then the shocks in this regime should 
be very weakly radiative. 

The obtained results are as follows. For internal shocks in e ± -pair-dominated 
outflows, the cooling time exceeds the microphysics time by many orders of magnitude, 
T CO oi 300 for the entire range of reasonable shock and outflow parameters. 
Therefore, emission from such pair-dominated shocks will have to be produced (if 
at all) in the far downstream region, not at the shock front. 

Somewhat similar conclusions follow for the external afterglow shocks in the 
electron-proton plasma. For most of the blast wave conditions, these shocks are also 
in the weakly radiative regime. The observed afterglow emission has to be produced 
in the large volume of the downstream region. The emission from the high-field region 
at the shock jump is strong for ultra-relativistic shocks with T > 350 in a high-density 
external medium, n cxt ~ 100 cm~ 3 , and for even higher T's at lower densities. Figure 
[4] shows Tcooi as a function of T. Dark and light blue regions indicate the range of T 
for which an external shock is in the radiative regime. 

The cooling time for the ISM and Wind models of the ambient density profile 
as a function of time after the burst is shown in figure [5] We obtained that, except 
for the very early times, the shocks (in both models) are characterized by very large 
T CO oi, indicating that the observed emission should be produced in the far downstream 
region. In the ISM model, the shock can be radiative only at the very early times, 
earlier than a second after the burst. In contrast, the shock in the Wind model can be 
in the radiative regime until about a hundred seconds after the explosion. Usually, the 
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afterglow sets in at least a few tens of seconds after the explosion. Therefore, there 
is a chance that the very early afterglow emission from a shock in the Wind medium 
is produced at the shock front itself. We speculate that one can expect a break in 
the afterglow (steepening of the light-curve) at about 10 to 100 seconds, indicating 
the transition from the radiatively efficient to the radiatively inefficient shock. It 
is likely, though, that this early afterglow emission will be swamped in a brighter 
prompt or high-latitude emission. It is far too early to make quantitative predictions, 
however, because radiative cooling can change the shock formation and evolution 
on a microscopic scale. A further study requires a radiative PIC simulation. It is 
quite possible that magnetic fields may be produced by the relativistic Richtmeyer- 



Meshkov instability (the vorticity-generating fluid instability, Sironi & Goodman 2007 
Goodman fc MacFadyen||2007 Milosavljevic et aT"||2007 ) in the afterglow phase. The 



dependence e B cx can be a benchmark of the model 

The most remarkable results are obtained for internal shocks in the baryon- 
dominated scenario. The "diagrams of state" - the contours of constant T coo i in 
the planes R-T and R-Ti are shown in figures [2] [3] One can see that the regime 
of strong cooling is limited to low bulk Lorentz factors, T < 60 and small radii 
R < 10 10 — 10 11 cm. Such shocks are well below the photosphere, located at 
i?ph ~ 10 14 — 10 15 cm for such low T's, which makes them hardly observable. However, 
radiative shocks with T coo i ~ 100 — 300 can form above the photosphere for the 
conventionally assumed GRB parameters, T ~ 200 — 300, L ~ 10 52 erg/s, Tj > 2.5, as 
is seen from figure [3] In this figure, we considered shocks forming at the photosphcric 
radius, i.e., when R4 = i?o(2r 2 ) ~ R P h- One should keep in mind that internal 
shocks can occur at radii smaller than i?j if the outflow modulation occurs due to 
instabilities in the outflow/jet itself, at radii R 3> Rq where the flow is already 
ultra-relativistic, rather than at the base of the jet, at R ~ Rq (Rees & Meszaros 



20051. This substantially relaxes the conditions on the optical depth, especially at 
high r. Note that depending on the efficiency of conversion of the kinetic energy 
into radiation, the radiation from internal shocks can produce the second, e ± -pair 
photosphere. Should this happen, the optical depth is estimated to be from few 
to ten. One can expect, therefore, the appearance of a thermal component in the 
spectrum. Such optical depths are, nevertheless, not enough to completely smear out 



a non-thermal component; t± greater than a hundred are likely needed for this (Pe'er 



& Waxman 20041. The pair opacity is less of a problem for low-energy GRBs: for 
L ~ 10 4S — 10 4y erg/s, T ~ 200 and the radiation efficiency of about 10%, the radiative 
shocks can occur at t± ~ 1 and well above the baryonic photosphere. 

The existence of the regime of GRBs with strong radiative cooling has many 
interesting implications. 

The first of them concerns radiation emitted at such shocks. The magnetic 
fields at the shock front are highly inhomogeneous and anisotropic on very small 
scales, much smaller than the typical electron gyro-radius. The fields predominantly 
have the filamentary structure, reminiscent of the filaments of the Weibel instability 



undergone subsequent numerous mergers in the foreshock region (see, e.g., Spitkovsky 



2008 Chang et al 2008). It has been predicted that radiation produced in such 



small-scale fields (called the "jitter radiation'' ) should be spectrally different from the 
conventional synchrotron radiation (Medvedev 2000[ ). The steeper than synchrotron 
low-energy spectral slope, F v oc v 1 , was an attractive solution to the "synchrotron line 
of death" problem in prompt GRBs ( Preece et al 1998 Medvedev 2000 1 . Because 



of the anisotropic, filamentary structure of magnetic fields at the shock front, jitter 
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radiation is predicted to be anisotropic too, with it's spectral shape varying with angle 
between the shock normal and the line of sight (Medvedev 2006). Combined with 



relativistic kinematics, the model reproduces certain spectral correlations observed in 
the time-resolved analysis of GRB d ata ( Kaneko et al |2006 1 , as well as the ubiquity of 



flat, F v oc v°, spectra in the sample (Pothapraga da et al |2007 1 . Afterglow lightcurves 



in the jitter regime have also been recently calculated (Medvedev et al 2007 Workman 



The second one deals with the fact that realistic GRB shocks in a certain 
parameter regime arc now accessible to simulations. Present 2D PIC simulations 
already resolve the shock formation and evolution, generation of magnetic fields in 
the foreshock region and their amplification toward the shock front, the electron 
acceleration to near equipartition with the protons and, tentatively, Fermi-type 
acceleration. The shocks in the radiative regime are remarkable in that the energetic 
electrons are able to radiate a substantial fraction of the shock energy, thus making the 
shocks radiatively efficient and, hence, observable as prompt GRBs. This opens the 
possibility to ultimately understand the properties of collisionless relativistic shocks. 
Moreover, by adding radiative cooling in the PIC simulations, one can study to what 
extent rapid cooling of electrons affects the shock structure on the microscopic scale. 
At last, one can obtain the radiation spectra directly from PIC simulations (this 
possibility has already been demonstrated by Hededal 20051, compare them with 



observational data and, after all, confirm or falsify the jitter model of prompt GRBs. 
All this does not seem feasible at present, however, for the shocks in the weakly 
radiative regime, in which the structure and dynamics of the magnetic field far behind 
the shock is crucial, but which is still difficult to trace with available computational 
resources. 

Third, in the radiatively efficient regime, an electron loses its energy every time it 
is crossing the shock, so Fermi acceleration of the electrons will be quite inefficient. One 
can expect that a power-law electron distribution will not form in this case and, as a 
consequence, the observed radiation spectrum may not have an extended non-thermal 
tail. Moreover, the peak of the emitted radiation is predicted for such a regime to be in 
the multi-MeV range, where the comptonized component is conventionally expected, 
but no strong second peak is expected in the keV-MeV range. This prediction can 
soon be tested with Fermi (former GLAST) observatory. 

Finally, the Weibel instability and radiation production in conditions very close 
to those in GRBs can be studied in existing laboratory experiments. The generation of 
filamentary magnetic fields, indicative of the Weibel instability, has been demonstrated 



in numerous laser-produced plasma experiments (see, e.g., Tatarakis et al 2003). In 
these experiments, a Petawatt laser beam with the intensity / ~ 10 2U — 10^ W cm - 
produces a relativistic plasma of density n ~ 10 19 — 10 21 cm -3 and relativistic electrons 
with the Lorentz factor of one to few hundred. Because the the density in the 
experiments is so close to the density in internal shocks, ~ 10 16 cm~ 3 , the plasma 
skin depths differ by about two orders of magnitude, cf., c/lu pp ~ 10 /im in the lab 
and ~ few mm in prompt GRBs. For the typical transverse size of a laser beam of 
~ 100 /j,m, a large number of the Weibel filaments are produced in the target, which 
makes the further study of the nonlinear dynamics of filaments feasible. By launching 
a probe relativistic electron beam through the target and performing the X-ray and 
beam diagnostics, as is proposed in the Hercules experimental] one can directly and 

1 http:/ /www. eecs.umich.edu/CUOS/research/index. html 
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simultaneously probe the magnetic field structure and the spectral properties of the 
emitted radiation. These experiments can be done in a number of existing laser 
facilities, such as Vulcan, Hercules, Omega, NIF and a few other. Apart from being 
a very exciting possibility for the Laboratory Plasma Astrophysics, such experimental 
studies are crucial for the verification of numerical results. 

In this paper, we evaluated the cooling time in GRBs with reasonable accuracy. 
One should understand, however, that there are intrinsic uncertainties in the values of 
certain parameters, which are difficult to quantify. In our calculation of the average 
residence time of an electron in a region of high magnetic field at the shock front, we 
did not take into account the effective filling factor of the field, nor did we consider 
the increase of the residence time due to the electron trapping. Such effects strongly 
depend on individual particle trajectories and are, therefore, sensitive to both the field 
structure and the particle distribution function. We, however, evaluated and included 
the effect of the dumpiness of the magnetic field on the effective (emission-weighted) 
£b, as it is relatively straightforward. Here we also used the results of 2D simulations, 
which may differ from a realistic 3D case. More importantly, in the strong cooling 
regime, the electron cooling occurs on the time-scale of the shock formation. This 
effect can, in principle, change the entire shock structure. Since the effect of radiative 
cooling has never been studied before, the results of this paper are suggestive. How 
does strong electron cooling change the collisionlcss shock formation, its structure and 
dynamics? The answer to this question can be obtained from future numerical PIC 
simulations which include radiative effects. However, we can speculate that shocks 
with strong electron cooling will likely have longer transition regions than the non- 
radiative shocks. This is because electron heating is an integral part of the electron- 
ion shocks ( Spitkovsky |2008 1 , as it enables electrons to escape the ion filaments. This 
reduces filament shielding, and facilitates ion filament merging, which leads to shock 
formation. If electrons are kept cold in the upstream, the shock may thus take longer 



to form (Lyubarsky & Eichler 2007 Gedalin et al 2008). Finally, in our study we 



completely ignored the magnetic fields present in the upstream region. These fields 
are relatively week, compared to the fields at the shock itself, but they can make 
a non-negligible contribution to the electron cooling if they occupy a large volume. 
Simulations of electron-positron shocks indicate that even in the longest simulations 
the upstream magnetic field has not reached a steady state, but continues to grow 



with the simulation time (Keshet et al 2008). In our analysis we also ignored the 



presence of the fields in the downstream, which have also been demonstrated to be 
affected by the upstream fields. Thus, our analysis is quite conservative since both the 
upstream and downstream fields can only decrease the electron cooling time and/or 
increase the radiative efficiency of the shock. 
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Figure 1. A snapshot of a shock from PIC simulation with T = 15 and 
mi/m e = 100. (a) and (b) — Maps of the particle density n and the normalized 
magnetic field strength eg, respectively. The incoming flow moves from right 
to left, while the shock propagates to the right. The simulation is performed 
in the downstream frame, (c) - Density profile averaged over the transverse (y) 
direction, (d) and (e) — The j/-averaged profiles of eg and e e (red curves) and 
the "emission-weighted" profiles of eg and e e (black curves) respectively. 
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Figure 2. Contours of T coo i vs. V for I\ = 4 (left panel) and a somewhat 
unrcalistically high T; = 10 (right panel) for the internal shocks in baryon- 
dominated outflows for T coo i = 10, 100, 300. Red filled regions correspond to 
^cool < 10 (the "strong cooling" regime), dark and light blue regions correspond 
to 10 < T coo i < 100 and 100 < T coo i < 300, respectively (the "radiative 
shock" regime), and the white region corresponds to the "weakly radiative shock" 
(Tcooi > 300). The black triangle outlines the range of parameters where radiative 
shocks occur and may be observable. In all cases, L52 = 1, L 7 = 0.1L, t v> —4 = 1. 
The radii, Ri, R p ^ and i? pa i r are plotted for reference. Internal shocks can occur 
above Ri line; optical depths, tj, and t±, are less than unity above lines il p h 
and -Rp a i r , respectively. White lines in the left panel denote radii at which the 
optical depths t), and t± are equal to 10 (dashed and dotted lines respectively). It 
is expected that some non-thermal spectral signatures of radiation emitted from 
not so large optical depths of t < 10 may still be observable in addition to the 
strong thermal component. 
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Figure 3. (left panel) Contours of T coo i vs. the Lorentz factor of internal shocks, 
Ti, for the internal shocks in the T = 300 ejecta. (right panel) Dimensionless 
cooling time, T coo i, in internal shocks at the baryonic photosphere vs. r\. Color 
coding and the outflow parameters are the same as in figure [2] 
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Figure 4. Cooling time in afterglows vs. the outflow Lorentz factor, for 
"cxt = 100 cm" . Color coding and the outflow parameters are the same as 
in figure [2] 
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Figure 5. Cooling time in afterglows vs. time after the burst, for the ISM 
and Wind models of the external density. In the left panel, we use the "typical" 
parameters E = 10 53 erg, A* = 1 and njsM = 1 cm -3 , whereas in the right 
panel, we use a rather extreme set of parameters, E = 10 54 erg, A, = 10 and 
n lSM = 100 cm" 3 . In both cases, we put a GRB at a typical z = 2; time in the 
plots ranges from 0.1 s to 100 days. Color coding is the same as in figure [2] 



